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ABSTRACT 


In  this  report  the  varianoe  of  narrow  band  noise  plua  c.w.  signal 
after  detection  and  gating  in  time  Is  derived.  Tho  problem  la  solved 
for  idoal  and  Gaussian  filters  with  square  law  and  linear  detection, 
the  latter  for  large  pvedeteetien  signal  to  noise  ratio.  In  addition  to 
those  filters  the  analysis  is  carriod  through  for  a model  which 
considers  the  envelope  of  narrow  band  noise  plus  signal  as  a function 
which  is  constant  for  intervals  of  duration  oqual  to  tho  reciprocal 
of  tho  noise  bandwidth  and  statistically  indopondont  in  different 
intervals.  A very  close  agroement  is  found  among  the  filter  types 
mentioned  and  tho  model.  Graphs  of  tho  mean  deviation  dividod  by  the 
mean  power  for  noise  alone  arc  given  for  both  filtors  and  tho  model, 
with  square  law  detection  and  predotoctien  signal  to  noise  ratios 
0,  0.5,  1,  2,  Finally  we  consider  a simplified  modol  which  is  similar 
to  the  model  proviously  mentioned  except  that  the  gate  is  assumed  to 
start  at  a point  of  change  in  the  model  voltage.  The  values  in  Tabic  1 
show  that  this  simplified  model  may  bo  used  **or  very  large  or  vory 
small  gate  widths  but  that  the  discrepancy  Is  considerable  for  gate 
widths  of  the  order  of  tho  reciprocal  of  tho  noi3o  bandwidth. 


Some  Statistical  Properties  of  Oated  Signal 
Plus  Narrow  Band  Noise 


1 


I.  The  Problem. 

This  paper  describes  some  of  the  statistical  properties  of  nar- 
row band  noise  plus  c.w,  signal  after  deteoiion  and  gating  in  time. 
More  specifically,  we  are  concerned  with  the  fluctuation  of  the  vari- 
ables f 

y(t)  » ( x( t)dt  (1) 


and 


t-6 


e(t) 


t 

j*  x(t)dt 
t-6 


t-T 

a 

x(t )dt 
t-T-6 


(2) 


(see  Figs,  la  and  lb) 

where  x(t)  la  the  output  of  the  system  shown  in  Fig.  2.* 


*(t) 


(C) 


(t  ) 


Fig.  1 


«The  noise  bandwidth  of  a filter  is  defined  as  usual  by 
where  Y(f)  is  the  transfer  characteristic  of  the  filter  . .>2 

and  fe  its  center  frequency.  Physically  j)  represents  lY  ' ro  • 1 

the  bandwidth  of  an  ideal  filter  which  has  the  same  total 
power  output  and  center  frequency  gain  as  the  given  filter. 


1 


2 


2 2 

The  variances  of  y and  s , which  we  denote  by  dy  and  c z , 

will  depend  upon  the  filter  and  detector  characteristics  and  are  in- 

2 

dependent  of  time,  since  we  shall  assume  x(t)  stationary,  dy  will 

2 

be  a function  of  the  parameters  p , P , w0  and  5 , and  d2  will 
in  addition  depend  on  the  parameter  T • 

We  shall  consider  the  following  band  pass  filter  characteristics! 

a)  Ideal  Filter  - transfer  characteristic  given  by 

Y(f)  * 1 f c - |-<f  <fc  + (3) 

a 0 elsewhere 

b)  Gaussian  Filter  - transfer  characteristic  given  by 

*( f-r0>2 

Y(f)  =.  e -~r-  (4) 

For  both  of  these  filters  0 <<  fc  and  the  noise  power  output  Is 

W - w0  J“  |Y(f ) |2df  - pw0 

In  addition  to  these  filters  we  consider  a model  for  the  envelope 
of  signal  plus  narrow  band  noise  which  is  frequently  used.  This  model 
considers  the  envelope  of  narrow  band  noise  plus  signal  as  a function 
which  is  constant  for  Intervals  of  duration  -4r  and  scatistioally  Jr.- 


3 


dependent  In  different  intervale.  An  example  of  euoh  a function  la 
shown  in  Fig.  3*  The  times  at  which  the  functions  changes  its  value 
ore  assumed  to  be  unknown. 


x(t) 


The  problem  is  solved  assuming  square  law  detection  and  for  large 
signal  to  noise  ratio  in  the  case  of  the  linear  detector. 

S.  0.  Rice  (Jour.  Acoustical  Soc.  Amer.  ll^  21 6 ( 19U-3 ) » referred 
to  nereafter  as  Ref.  A,  has  solved  this  problem  in  essence  for  the 
case  of  zero  signal.  The  present  work  extends  that  of  Rice  to  in- 
clude signal  as  woll  as  noise.  The  approach  is,  however,  different 
from  that  of  Rice, 

II • Noise  in  Linear  Systoms. 

The  gating  process  is  a linear  one  and  therefore  tho  theory  of 
noise  in  linear  systems  may  be  applied. 

If  we  consider  a linear  system  defined  by  a transfer  character- 
istic Y(f)  or  equivalently  by  a weighting  function# 


#Tiie  weighting  function  is  the  response  to  a unit  impulse  excitation. 


k 


f 2%  ift 
h(t)  » Y(f)e  df 


then  the  output  g2(t)  is  related  to  the  input  £^(t)  by  tho  convo- 
lution integral 


g2(t)  a j h(  T)gl(t-  X )dx 


s . that 


e22(t) 


■LI 


h(  x )h(  \)gi(t-  T)g1(t-  X)dx  dX 


g (t)  = i h(  x )h(  X)  9j(x  - X)dx  d\ 


v/htre  9^(x)  i3  the  autocorrelation  function  associated  with  g1(t). 
For  tho  gating  process  described  by  Eqn.  (1)  and  (2)  tho  corres- 
ponding weighting  functions  are 


hy(t) 


1 0 <t  < b 

0 elsewhere 


1 0<t<6 

h2(t)  a -1  T < t <T  +6  (7) 

0 elsewhere 

The  doublo  integral  in  Eq,  (£)  may  be  reduced  to  a single  integral  by 
the  transformation  x - X = t)  f X s { , 

For  h_(t)  as  given  in  Eq,  (6),  we  then  havo 

i rb  *mrl 

h(x)h(X)  9x(x-X)dx  dX  =2\  dn^H)  h(  f|  + Oh(  £)dS 

Jo  Jo  J0  Jo 

From  Fig,  If  for  hy<T|  + 0 and  hy(  K ) wc  have 

h(n  + Oh(  K ><u  * 6-n 


o 


for  0 a r)  a 6 


Pig»  k 

kio  that  the  double  Integral  of  Eq.  (5)  becomes 


2 r 

&2  (t)  = 2 (6  - q)  9X(  q )dq 


(see  footnote  to  Eq.  (9)) 

For  hz(t)  as  given  in  Eq.  (7)  wo  have 


,T+6  _T+  5 


o Jo 


h(  x)h(  X)^(T-X)dTd\ 


I 

^(q)  h(q  + 0 h(0  d ; 

n> 


From  Fig.  5 for  hz(q+5)  and  hz(0  one  may  note  that  there  are  the 

-T+6-q 

following  contributions  to  the  integral  h(q+£)h(Od£  : 

Jo 

a)  For  0 £ q = 6 one  has 

6-  q + (T+6-  q)  - T 3 2(  6 - q) 

b)  For  0 = T -q  = (i.e.  T - 6 = q = T ) one  has 

- [5  - (T  - q)  ] = [q  - (T  - 6)] 

o)  For  0 = T + 6 - q = 5 (i.e.  T a q s T + 6)  one  has 

- [T  ♦ 6 - q] 
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Pig.  5 


Thus  the  double  Integral  of  Eq.  (5)  becomes,  for  the  doublo  gate; 

b rT 

[q-(T-b)]  9i(n)dl 


g22(t)  * (6-n)<Pi(|J)<in  • 2 C«i 

i U 

T+6 

- 2 J (T+b-  t))  9j(  n)dn 


(9) 


This  result  Is  not  altered  In  the  event  that  the  regions  mentioned  in 

a)  and  b)  overlap,  _____ 

2 

It  may  bo  noted  that  wo  havo  chosen  to  express  g2  (t)  as  a 

doublo  integral  involving  h(t)  and  qp(  t)  which  is  then  reduced  to 

— 

a single  integral.  Wo  could  equally  well  have  expressed  g2  (t)  di- 
roctly  as  a single  Integral  involving  the  transfer  function  Y(f)  and. 
the  input  spectral  density.  This  procedure  was  not  followed,  howevor, 
since  the  integrations  to  be  performed  aro  not  nearly  as  simple  as 

*VThe  variance  <s  ^ s gg^(t)  - [g2(t)]  ia  givon  by  Eqs,  (8)  and  (9). 
with  replaced  by  minus  its  d.c,  voluo. 


uO 


7 


those  indicated  In  (8)  and  (9)  for  tho  case  of  tho  Oausslan  filter 

and  the  sxxlel,  We  thus  uso  (8)  and  (9)  for  the  evaluation  of 
2 

and  *s  respectively  (soe  footnote  to  (8)  and  (9))»  for  whioh  we 
need  the  autocorrelation  function  associated  with  the  output  of  the 
detector,  This  Is  discussed  In  the  following  section* 


III*  Properties  of  the  Detector  Output* 

The  action  of  the  detector  1s  such  that  It  ratifies  linearly 
the  Input  signal  plus  narrow  band  noise,  raises  the  result  to  some 
power  and  then  filters  so  that  all  components  of  frequency  of 
the  order  of  tho  carrier  frequency  and  above  are  stopped*  In  our 
analysis,  however,  we  consider  the  detector  as  a device  which  first 
obtains  tho  envelope  of  the  input  signal  plus  narrow  band  noise  and 
then  raises  tho  envelope  to  tho  power.  This  output  is  the  same 
as  that  of  the  actual  detector  apart  from  a constant  multiplying  fac- 
gor,  (See  Pig,  6 below,)  For  purposos  of  computation  It  is  more 
convenient  to  consider  the  second  description  of  detector  action. 


signal  plus 
narrow  band 
noise 


signal  plus 
narrow  band 
noise 


x(t) 


x(t) 


Pig*  6, 
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The  speotrum  and  the  equivalent  autocorrelation  funotion  of  tho 
output  of  a square  law  dotootor  have  boen  obtained  by  Rice  (Boll 
System  Tech,  Jour.  2j^,  Part  IV  (1945))  and  by  Lawson  and  TJhlenbeck 
{ Threshold  Signals,  Vol,  2lj.  of  M.I.T.  Rad.  Lab.  Ser.,  McGraw-Hill, 
1(;50.  Sec.  7,2,  referred  to  hereafter  as  Ref,  B),  Tho  result  of  In- 
terest for  the  determination  of  0 and  0 2 is  the  autocorrola- 

J * 

function  minus  its  d.c.  value.  (See  footnote  following  Eqs* 

'")  and  ^9)).  This  is  given  in  Ref.  B (p.  155,  Eq.  (17))  as 

q>(  t)  = IjA  p(  x)  + 4w2  p2(  x)  (io) 


W = noise  power  input  = fiw0 
P = signal  amplitude 


and 


p(  t)  * 


1 | Y | 2 cos  2 w(f-fc) x df 
J o 

[ I y|  2 df 

Jo 


(11) 


In  the  case  of  linear  detection  tho  autocorrelation  function  may 
rlsa  bo  obtained  from  Ref.  B (p.  1$7,  Eqs.  (20a)  nnd(20b)),  but  in 
toe  case  of  signal  much  groator  than  noise  it  may  bo  obtainod  roadily 
by  examining  the  Rico  representation  of  narrow  band  noiso  as 

X(t)  sin  2 w fct  + Y(t)  cos  2 wf0t 

wbox-o  X(t)  and  Y(t)  are  slowly  varying  functions  which  have  Gaussian 
probability  distributions  with  variance  W and  have  normalised  auto- 
correlation functions  = p(  x).  Thus  tho  envelope  of  narrow  band 
nriso  plus  signal  * p cos  2wfct  may  be  represented  at  Rny  instant 
by  the  vector  diagram  shown  in  Pig,  7,  and  for  P2>>x2  ■ Y?  the  fl  v • 
tn a'; ion  in  R is  duo  primarily  to  tho  fluctuation  in  X so  that  in  tho 
Unit  of  very  largo  P,  the  autocorrelation  function  minus  its  d.c, 
value  is 
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R(t)  R(t  + t)  -(W)2  » X(t)  X(t  + x)  rn  pw 
Thus,  for  the  Ideal  filter,  substituting  (3)  in  UD* 


J 


p(  t)  - V0/2 


eoa  2*  (f-fc)  Tdf 


df 

Vf>/2 


(12) 


sin  n t B 

*xp 


(13) 


ard  for  the  Gaussian  filter,  substituting  (lj ) in  (11) 

i(f-fc)2 

2 cos  2n(f-f  ) x df 
P c 


p ( t ) = 


,OQ  P 

- * (f-fc) 
o 3 — df 


T 


,O0 


- oo 


- nf* 

o cos  2 nf  xdf  _ „ 

P 02  2 

- It  0 T 

= e C • ) 


r.- 


Hf2 

l7 


in -jo  -2.  >>1.  Eqs.  (13)  and  04)  may  thon  bo  substituted  in  (10’  to 

P 


PlG*  7 
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obtain  the  autocorrelation  function  mlnua  its  d.c,  value  in  the  oaae 
of  the  Ideal  and  Gaussian  filters. 

For  the  model  the  autocorrelation  function  is  given  in  Fig*  8 
in  the  case  of  the  square  law  deteotor. 


Fig,  8 


From  Ref,  B (p,  155,  Eqs,  (lljb)  and  (17))  we  have 
2 

X^(t)  - X2(t)  » 4WP2  ♦ IfT/2 


so  that  the  autocorrelation  function  minus  its  d.c,  value  is  given  by 

(15) 


*(  t)  » U^WP2  + ^)(1  - p K|  ) |T|  <A 

P 


M >“jr 


In  the  case  of  the  linear  doteetor  the  autocorrelation  function  is 
given  in  Fig,  9*  * 2 

For  P^»W,  x^(tl  - x(t)  may  bo  obtained  either  from  Ref,  B 
(p,  155,  Eqs,  (ll|A)  and  (lifb),  retaining  the  first  two  terms  in  the 
asymptotio  expansion  of  (ilia))  or  directly  by  sotting  x » 0 in  (12), 
giving  2 

ViT)  - x(t)  ■ w 


11 


X Tt) 


so  that  the  autocorrelation  function  minus  itr  d.c,  value  is 

9 ( i)  ■ W(1  - 0 hi) 


hi  S I 


(16) 


= o 


hi  > \ 


IV*  Results,  Discussion  and  Conclusions, 

2 2 

<j,r  and  o can  now  be  determined  for  the  ideal  and  Gaus^Jan 

y z 

filters  from  (13),  04),  (10)  and  (12)  by  substituting  in  (8)  and  (9), 
and  for  the  model  by  substituting  05)  and  06)  in  (8)  and  (9),  The 
integrations  are  carried  out  in  the  appendix  and  the  results  are  as 
follows : 

For  the  single  gate,  square  law  detector 

.2  f0(2  * p 6) 


Ideal  filter  0, 


2 __ 

y s ( * p )2 


rG(2  * p&)  . , A*  , % 

j +4nwP&0'(*  P 6)J  07) 


•MTo  perform  the  Integrations  it  was  found  convenient  to  introduce  the 
function  x * , 

0(x)  -f  ds  £ dt 


defined  by 


on(*)  - — -ff-2- 


0*(o)  * G(o)  = 0 * 


< < 

The  values  of  G(x)  and  G'(x)  for  0 ■ x * 100  are  given  in  Table  2 at 
the  end  of  this  report* 
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Gaussian  filter 
2 W2 


2- Tsfift  SZn  Erf(06  ]/Fn  ) + -~r  e 

W3T  0 2 \/k  L V* 

Model  d 2 = IfW2  52(1  + 2 |i)(l  - -££-),  06  ^ 1 


2 .-2  *0262  _ 2 


1, 


= l*W2  62(1  + 2 *i)  .f  ? - " ■!■■  , 06  = 1 


(18) 

(19) 


(06)‘ 
p2 

wliore-  ii  = i-  * predetection  signal  to  noise  ratio, 

2W 

2 

For  06  «1,  0y  for  both  filters  as  well  as  the  model  should 
be  the  variance  of  a single  sample,  of  width  6 , i,e»  we  should  ex- 
pi  -it.  to  have 


c^(t)  - x2(t) 


= 6 2[ ifV/P2  + 4w2  ]=  IfV/2  62(1  + 2 n) 


(20) 


(see  Ref.  B.,  p,  1*?£,  Eqs.  (Tij-b)  and  (17)).  Substituting  the  expres- 
sions for  G(x)  and  Erf(x)  as  given  in  the  appendix  (Eqs.  (A-ll)),  (17) - 
(Ifi)  and  (19)  reduce  to  (20)  in  the  limit  for  06  « 1, 

For  06>>1»  we  mav  use  the  asymptotic  expressions  for  G(x)  mid 
Err (x)  given  in  the  appendix,  Eqs,  (A-12),  and  obtain  the  following 
limits  from  (17,  ( 16 ) and  (19): 


Ideal  filter 

\ - *5  u 
7 * 

♦ 2 ti)06 

(21) 

Gaussian  filter 

»% 

f + 2 jiy  06 

(22) 

y 0 c \ * 

7 

Model 

- *5 

7 0* 

♦ 2 |i)  06 

(23) 

In  each  of  these  cases  dy  is  proportional  to  6 as  is  to  be  expected 
since  for  06  » 1 the  gating  process  is  effectively  adding  up  a 
nvT.oer  (of  order  6 ) of  independent  random  variables,  the  varianoe  of 
the  sum  then  being  the  sum  of  the  variances.  It  is  to  be  noted  that 


U 


both  for  (36  « 1 and  (36  >>1  the  0 corresponding  to  the  ideal 

filter  and  the  model  are  equal  for  all  signal  to  noise  ratios*  The 
p 

values  of  d y for  the  ideal  filter  and  model  do  not  differ  appreci- 
ably even  for  (36  *«•  1,  as  is  shown  by  the  figures  given  in  Table  1« 

For  the  particular  case  of  zero  signal  the  expressions  given  in 
(17)  and  (18)  for  d y2  for  the  ideal  and  Gaussian  filters  with  square 
law  detector  and  single  gate  are  compared  with  the  corresponding  ex- 
pressions given  in  Ref.  A (p*  22 3,  Eq*  (6,10)  and  p,  225>)  in  the  ap- 
pendix, and  are  shown  to  be  four  times  the  variancos  given  in  Ref.  A* 
Th.<  s is  to  be  expected  since  in  Ref.  A Rice  considers  the  short  time 
average  value  of  the  noise  after  detection,  whereas  we  have  been  con- 
sidering the  envelope,  and  those  differ  by  a factor  of  two# 

For  the  single  gate  and  linear  detector  with  largo  signal  to 
no? 3e  ratio,  we  have  the  following  results: 

Ideal  filter  dy2  = 2W  6 ) 

Gaussian  filter 
0 2,— rA 

y 2p  ‘w* 

Model  dy2  =s  W 


Again,  for  (3  6 < < 1,  d " for  both  filters  as  well  as  the  model 
should  be  the  variance  of  a single  sample,  of  width  6 , i.o#  we 
should  expect  to  have 

dy2  * 62[x2(t)  - xUl]**2*  (27) 
Substituting  the  expressions  for  G(x)  and  Erf(x)  as  given  in  the  ap- 


f / o -*P2&2  J'\ 

|2  (36 /k  Erf  (p  bjn  ) 0 " "v'T  J 


>2(1 

■ - 1/3) 

( pb  )2 
2 


p6=  1 


p6  = 1 


(2  b) 
(26) 


pendix  (Eqs.  (A-ll)),  (2i<.),  (25)  and  (26)  each  reduce  to  (2?)  in  the 
Unit  for  05  «1, 

For  05  > >1,  we  may  use  the  asymptotic  expressions  for  Q(x)  and 
Erf(x)  given  in  the  appendix,  Eqs.  *12).  We  then  obtain,  from  (2!;), 
(25)  and  (26), 

aj  = • £6  (27») 

* r 

f^i'  both  the  ideal  and  Gaussian  filters  and  the  model.  In  (2’jh),  as 
in  (21),  (22)  and  (23),  <Jy2  is  proportional  to  5 , as  is  to  be  ex- 
pected. 

For  the  double  gate  with  square  law  detector  wo  have  tho  follow- 
ins-  results: 

Ideal  filter 

G(2  x 0 5)  + 1+Ux0  & G»(  x 06  ) * 

+G(2  x 0 T)  + 1*.^  x 0 T G'(*  0T) 

(2  r. 

G(2  x 0 ( T—  5))  + lj.|ix0  (T-  5)G'(x9  (T-6))] 

l 

I 

+G(2  x 0(T+  5 ) ) + x 0 (T+  5 )G»(  x0  (T+  5 : j 


Gaussian  filter 
^ 2 2W2 


02/x 


2 0 5 Erf (05  VcTx  ) + 20  T \/2x  Erf(0  T /Zr,  ) 

-0  (T-  5 ) >/2xErf(  0 (T-  5)  y£x  ) - 0 

- 0(T+  5)  ■/Z*Ert(  0(T+  5)  \£x  ) 

2 -2  x02  52  2 

+ /5t  6 . " 

-2x  0V 


1 f-i-  “2^2(T-  „_2_ 

2 L*  o ;j-  vn 


2 -2x02(T+5)2 

+ e 

Vx 


(29) 


> 1 
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Gauaaian  filter  (cont*d) 

2p6  Erf  ( p 6 /*  ) + 2 p T Erf  ( p T /%  ) 

QjLBL-,  1 - P(T-  &)>/*  Erf(P  (T-  6 )v/Jc  ) - p(T+i  y«Epf(  p (T+i^/if 


P2*'* 


2 - n P 6 


2x  2 


+ 7T  ° 


JL 

/? 


! f 2 -*P2(T-6  ) 

-il^T0 


u oT 

2 - a P2(T+  6 )2 

♦ -£r-  e J 

v/n 


2 - *P2T2 

■2C^° 


Model 
a)  d ' 2 = 


8VV2  (1+2  jt  )(  p&  )2 

P7 


(JT-i 


P (T+  6 ) < 1 


b)  d 


2 


(1+2  p ) 

r 


( Pb  )2  - P • B2(T-  b )*~  + 0^(T-  l 

2 e~LT 


- + 1 + ( 0 T)2  - I(  p T)3 


o)  d 


8w2 

Is 


pT  <1  - 0 (T+  6 ) 

f< ?» >2  .^»>)3.i.eittsa!*ei4=jg 

(X+2|i  ) L J 6 2 2 


(22.) 

3 


P 


(T-  6 ) 


P&  < 1 

P(T-  6 ) < 1 i p T 


8W2 


d)  d * » ~j(l+2^  ) [ ( pb  )2  - j(  p 6 )3  ] 


P 6 < 1 

1 = p (T- 5) 


' ''  *0*  **  ■ *“**'  - 
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«)  «,a-^s<uatl) 


PlZr_fi2  + s2(t- 


n o. 


8W2 


£&+2|i)(P  b - i) 


lip* 

P (T-  b)  < 1 S PT 
1 * p6 

1 i P (T-  b) 


For  P&  « 1,  dyW  should  be  the  variance  of  the  difference  of  two 
samples,  each  of  width  6,  separated  by  a distance  T,  i,e.  we  should 
erpect  to  have 

. 2 _ 


(x2(t)  - 


x2(t+T))2  - (x2(t)  - x2(t'+  T ) j 


a 2 6 2 [ x^-(  t ) - x2(  t )x2  ( t + T ) ] 

a 8w 2&2  C(l-p  2{T) ) ♦ 2jl  (1-  p (T) ) 3 
(See  Ref.  B,  p.  155,  Eq.  (17)) 

For  the  ideal  filter  this  is 


(31) 


2 2 2 
0 m a ew  bd 

Z 


Fi  . zlrts-U  + z»(i.  sis.  -al 
L ( x P T)2  ' xP  T / 


a > 


and  for  the  Gaussian  filter 

o!2  = ew2»2(i.e-2"P2l2t2ll(i.,-"f|2ia)] 


(3ib) 


(31a)  and  (31b)  may  be  obtained  directly  from  (28)  and  (29)  by  examin- 
ing these  expressions  in  the  limit  of  P6«l, 

In  the  case  of  the  model  we  have,  for  P&«1 


8W252(1 

+ 2 H ) p T 

for  pT  < 1, 

from  (30a) 

(31c) 

8W2&2(1 

+ 2 M- ) 

for  PT  > 1, 

from  (30d) 

(3H- 

For  p&  «1,  PT  >>  1 we  have,  for  both  filters  as  well  as  fjr 


*eifc»v 
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the  model, 

d82  « Qw2  62(1  + 2 *0  (32) 

which  is  twice  the  variance  of  a single  sample  of  width  6 (see  Eq. 
20),  as  is  to  be  expected# 

For  p&  » 1,  we  have,  substituting  the  asymptotic  expressions 
for  G(x),  G’(x)  and  Erf(x)  in  (28)  and  (29),  and  directly  from  (3°)» 

« 2 = ^4(1  ♦ 2(i)P&  (33) 

* P2 

for  the  ideal  filter  and  the  model,  and 

'A) 

2 

0 as  given  in  (33)  and  (3^)  are  twice  the  variance  for  a 
z 

single  gate  for  which  P*>  >>  1 (see  Eqs,  (21),  (22)  and  (23)),  as 
should  be  expected.  It  may  be  noted  that  (33)  and  (3t)  are  obtained 
for  P6  » 1 regardless  of  whether  p(T-  6)  « 1,  p(T-6  ) 1 or 

P (T-  5 ) » 1. 

For  the  particular  case  of  zero  signal  the  expression  given  in 
(27)  for  dy2  for  the  ideal  filter  with  square  law  detector  and 
double  gate  is  compared  with  the  corresponding  expression  given  in 
fief.  A (-«  22k,  Eq.  (7.3))  in  the  appendix,  and  is  four  times  the 
variance  given  in  Ref.  A,  as  it  should  be,  since  in  Ref,  A the  aver- 
age value  of  the  noise  rather  than  its  envolope  is  considered. 

For  the  double  gate  and  linear  detector  with  large  signal  to 
noise  ratio,  the  results  are  as  follows: 

Ideal  filter: 

« 2 . -M—  [*B60'(  *P»  ) + * P T0'(  *P  T)] 

1 Up)2  (3  S) 

. w 

<*P>2 


[xp  (T- 6)Q’(x  P (T-  b))  + xp(T+& )0»(*p(T+  f* 


08 


Gaussian  filter * 

m 2 W 


P2/^ 


2 p6  t/n  Erf(  p&|/5  ) ♦ 2 PTi^T  Erf(p  T l/x  ) 

- P (T-  b )l/5  Erf  ( p (T-  6 )i/%  ) - 

- p(T*  b)SH  Erf(  P (T+6  )/*  ) 


2 -*P252  _£_ 

vie  ® ’ /* 


- 1 F 2 e“*P2(T- 

2 Lj/5 


6 )2  . 


2 -np2T2  2 -*P2(T+ & )2I 

-2  7Te  +7Te  J 


(36) 


Model: 


a)  dj2  = ^■"’5,(  P&)2  [pT  - yPb]  p(T+  6)  < 1 

P 


b)  <J 


2W 

7* 


( P6  )2  - i(  P6  >3  - pZ(T-  fc)2  + B3(T-  &)3 
3 2 6 

- pi|±_£l  + I + ( pT)2  - i(  pT)3 


PT  < 1 - P (T+6) 


% 2 21V 

o)  d«  -72 


( P6  )2  - i<  ft  )3  - i . p2(»-  6 )2  P3(T-  6 >3; 


lf-v, 


♦ B(T-  6) 
2 


d)  d 


2 = -^§  C(  81  )2  - i(  P&)3] 

a P<*  3 


p5  < 1 

P(T-  6 ) <1  i pT 

p5  < 1 

1 = P (T-  6 ) 


e)  0, 


-fjp- 1 - c2^  ♦ 


p6 


P ( T-  6)  < 1 i pT 


t ) c.‘  - -f?<  P6-  i) 


2 2W, 

8S  " 

■ P 


1 * P& 

1 = p(T-  6) 
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Again,  for  (3  b <<1,  <j  2 ahould  be  the  variance  of  tho  differ.* 
ence  of  two  samples,  each  of  width  separated  by  a distance  T, 
i.o.  we  should  oxpect  to  have  2 

d22  = 62  [(x(t)  - x(t+T))*  - (x(t)  - ”x(t+T))l 

« 2b  2 [x2(t)  - x(t)x(t+T)] 

« 2 bhfil  - p(T)  )#  # OB) 

f^r  large  signal  to  noise  ratio* 

For  tho  ideal  filter  this  is 

0 2 = 2W6  2(1  - 2l2L*il)  (38a) 

z *p  T 

and  for  the  Gaussian  filter 

HzZ  > 2V!t>  2(1  - e"*^  2t2j  (38b) 

(38a)  and  (38b)  may  be  obtained  directly  from  (35)  and  (36)  by  exanir. 
ing  these  expressions  in  tho  limit  of  Pb  < < 1* 

In  the  case  of  the  model,  we  have,  for  Pb  < < 1, 

d2  = 2Wb2«pT  for  p T < 1 , from  (37a)  (30; 

z 

0 Z » 2W  b 2 for  pT  > 1 , from  (37d) 

For  pb  » , pT  > >1,  we  have,  for  both  filters  as  well  as  for  the 
model 

•»From  Eq.  (12),  for  large  signal  to  noise  ratio,  wo  have  2 

x(t)x(t+T)  - x(t)  « p (T)W,  and,  setting  T ■ 0,  x2(t)  - x(t)  ■ W. 

Thus  x^(t)  - x(t)x(t+T)  » W(l-  p(T) ) 
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e2«2IT&2  (39) 

whieh  la  twlee  the  varlanee  of  a single  sample  of  width  b (see 
Eq.  (27). 

# 

For  j)b  > > 1 we  have,  substituting  the  asymptotic  expressions 

for  G(x),  0*(x)  and  Erf(x)  in  (3$)  and  (36)  and  directly  from  (37), 


(1*0) 


a 2 .2 
* P 

for  both  filters  as  well  as  for  the  model.  The  remarks  made  follow- 
ing (3^)  apply  equally  well  to  (1+0),  (of.  Eq.  (27b). 

In  the  accompanying  graphs  (Figs,  10,  11,  12),  the  following 
plots  have  been  drawn,  in  each  case  for  the  square  law  detector  only: 

v dv 

1)  -=*  as  a function  of  06  for  both  filters  and  the  model,  and 

p2 

in  each  case  for  pre-dotection  signal  to  noise  ratio  p = — = 0, 

cW 

0 5,  1,0  and  2.0,  v/here  m = 2W  6 in  the  mean  power  for  noise  alone 
and  gate  width  6 . This  same  quantity  has  been  plotted  in  Ref.  A 
(p,  217,  Fig.  1)  for  the  ideal  and  Gaussian  filters  in  the  case  of 

C m 

n**  signal,  and  is  denoted  there  by  i • 

„ • "t 

2)  — 5.  as  a function  of  B5  , with  T * b , for  both  filters  an.' 

in  o 

r 

uhe  model,  and  in  each  case  for  signal  to  noise  ratio  p = — » 0, 
0,5 » 1*0  and  2.0,  whore  m *»  2W  6 . For  the  particular  case  of  the 
ideal  filter  and  no  signal  this  plot  is  givon  in  Ref.  A (p.  217,  Tig, 


1)  and  is  denoted  there  by 


JT,T 


nvp 


3)  as  a function  of  0 T,  in  the  limit  as  6 — **  0,  for  both 
m r • 

f 5 Iters  and  the  model,,  with  m and  tho  values  of  +1  as  given  in  above. 
These  graphs  show  a very  good  agreement  among  the  filter  types  chosen 
find  the  model.  For  most  applications  use  of  the  model  is  probably 
justified.  It  should  be  noted,  however,  that  when  the  model  is  used 
in  a very  simple  manner  for  computations,  that  ie,  by  assuming  that 


21 


the  gate  starts  at  a point  of  ohange  In  the  model  voltage,  then  the 

results  are  not  as  good,  as  is  shown  by  the  values  of  J!z  given  In 

* m 

Table  1 for  06  * 0,  1,  2,  3,  Ij.,  *>• 

Applying  the  model  to  the  gating  process  is  essentially  the  saire 
as  adding  a oertaln  number  of  Independent  variables*  When  the  model 
la  used  in  the  simple  manner  desoribed  above,  the  number  of  inde- 
pendent variables  added  is  neoessarlly  less  than  or  equal  to  the  num- 
ocr  added  when  the  model  is  used  as  in  this  paper*  For  equal  gate 
durations,  both  methods  will  give  the  same  mean  value,  but  the  simple 
method  will  result  in  a larger  variance,*  For  06  > > 1 and  for 
0*>  < < 1,  both  methods  agree  and  the  greatest  difference  Is  whor*. 

06  ^ 1. 


>;-As  a simple  example  of  this  principle,  consider  two  Independent  i on- 
don  variables  x and  y having  the  same  variance  oz.  If,  however, 
compute  the  variance  of  x+v  we  get  o2  , 


DOUBLE  GATE  Hi  vs;  /3 £~:um:.TA&z 


■fi  i 


I f * fiS 


i i rt"*0 


o 

i .;T- 


i::;u 
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TABUt  I 

Values  of  4t  for  p ■ 0,  0,5#  lt  2 


P&  * 

0 

06  « 1 

fib  m 2 

fib  - 3 

p&*  4 

fib»  $ 

; _ i*-0 

Simplified  Model 

1 

1 

0.5000 

0.3333 

0.2500 

0.2000 

tfudel 

1 

0.6667 

0.4167 

0.2963 

0.2292 

O.I867 

1 Ideal 

1 

0.6559 

0.3961 

0.2825 

0.2196 

0.1797 

1 

i Oausaian 

1 

0.5479 

0.3225 

0.2184 

0.1665 

0.1341 

i ■ 

I o*5 

f 

i Simplified  Model 

2 

2 

1.0000 

0.6667 

0.5000 

0.4000 

1 M'del 
1 

2 

1.3333 

0.8333 

0.5926 

0.4583 

0.3733 

| Ideal 

2 

1.1*296 

0.8476 

0.5929 

0.4571 

0.3716 

i 

i l.auaaian 

2 

1.2312 

0.7431 

0.5165 

0.3968 

0.3217 

H * 1.0 

i Simplified  Model 

3 

3 

1.0000 

0.7500 

0.6000  j 

| i iodel 

I 

3 

2.0000 

1.2500 

0,8889 

0.6875 

0.5600 

j 'deal 

3 

2.2033 

1.2990 

O.9032 

0.6946 

0.5634 

1 rauaaian 

3 

1.911+4 

1.1637 

0.8145 

0.6272 

0.5092  | 

i 

L a.o 1 

! Simplified  Modil  5 

1 

5 

2.5000 

1,6667 

1.2500 

1.0000 

Model 

5 

3.3333 

2. 0833 

1.4815 

1.1458 

0.9333  i 

Tdaal 

5 

3.7507 

2.2019 

1.5239 

1.1696 

0.9471 

Gauaalan 

5 

3.2809 

2.0049 

1.4106 

1.0879 

0.8843  | 
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In  this  appendix  we  give  the  detalle  of  the  integrations  in- 

2 2 

volved  in  the  determination  of  d and  d , 

Let  ue  consider  first  the  ideal  filter*  Substituting  (10)  and 
(13)  in  (8)  we  have,  for  the  square  law  deteotor,  single  gate. 


= 2^* (6  - t ) p, 


,2  sin  x 8 x 
X P T 


* (^JJ 


(A-l) 


To  perform  the  integrations  it  was  found  convenient  to  Introduce  the 


function 


G(x)  * J ds  j — p-  — 


(A-2? 


defined  by 


o»<*>  = issp 

X* 


G'(O)  » 0(o)  = 0 


(A-3> 


and  in  terms  of  which 


SiS-S  a 2G"(x)  ♦ xG"(x) 
sin2  | 2(1  . oos  x) 


(it 


2G"(x)  . 


Thus,  nuking  the  necessary  substitution  of  variables  in  (A-l)  we  have 
2 *>6 

<J  2 a §££-  j (ups  »x ) [ 20*  ( x ) + xG*  ( x ) 3 dx  (A-.’p 

7 (*p)2  / 


42f_  ft 

«*  >2  J - 


xpb 

(2  *06  - x)Gn(x)ds 


Integrating  successively  by  parts  and  noting  0*(o)  a G(o)  * 0 we  then 
obtain 


y ( x^ )‘ 


[0(2  X06  ) + 4 o»(  * M ) 1 


U-S) 


where  u * — — • Is  the  pre-detection  signal  to  noise  power  ratio*  The 
2w 

result  for  the  ideal  filter  with  linear  detector,  single  gate  and 
U > > 1 is  obtained  similarly,  substituting  (12)  and  (13)  In  (8), 
giving 


■ 2 ! (&  - 


x )W 


sin  npx 


o 

2W 


KpT 

wpb  0 • ( x p b ) 


dx 


(a-6: 


( * p )2 

1516  results  for  the  ideal  filter  with  the  double  gate  may 

bo  obtained  in  similar  fashion*  Substituting  (10)  and  (13)  in  (9), 

• * 

we  have,  for  the  square  law  detector, 

= *&£■  * #*  ( d, 

o 

h-  (t-  b )]  [4WP2  ♦ 4W2  (g&ffi-  )2]  « (A--,) 


* w 

-2  j [T+b  - xK4*p2  — ■ ■ ♦ 4W2  (sAYpV  ) 3 


d^ 


T 

8W2 


0(2  xpx  ) + 1^*  p b 0*(  it  P 6 ) 
(x  P)*  |_+  0(2  T)  + ^ xP  TO  * ( it  P T ) 

liW2 


0(2  tt  p (T-  6))  + Ij.uk  P (T-  &)G«(  x P (T-  6)) 

( xp  )c  |+  G(2  up  (T+  fc)  ♦ Ij-U*  P (T+&  )G»(x  P (T+  6)) 
For  the  linear  deteotor,  ^ » 1,  we  have,  substituting  (12)  and 
(13)  in  (9) 


*,2  -*[*- 


t]  waiJiaSi-  4, 

xfix 


-2  /[,-(►*« 


d * 


T-6 

T* 


-2  [[**-, ].**?*-*, 


2$ 


(A-S) 


a -!*£ [wP60»(  xpb  ) + *p  T G*(  up  T)1 

(up  )2 


2W 


UP  (T-  6)G’UP  (T-  6))  xfi  (T+  &)G»(*P  (T+  6)3 


UP  )2 

We  consider  next  the  Gaussian  filter*  Substituting  (10)  and 
(l)j.)  In  (8)  we  have,  for  the  square  law  detector,  single  gate 

2 ( 6 • t)  ftWF2  A .-2*e2*  2]  4, 

(A-9) 


(8)  we 


B6/x  2 B6V57 

^~2  J S*  -x)e’x  dx  + f (fib  '/2x  - x)e 


-*2 


dx 


0 o 

where  we  have  made  the  substitution  x * fix  l/x  to  obtain  the  first 
Integral  and  x * Px  \/2x  to  obtain  the  second.  The  Integrations 
are  now  straightforward  and  may  be  expressed  in  terms  of  the  funs ‘Jen 

.2 


Brf(x)  « 


dt 


(A-iO) 


The  results  for  the  Gaussian  filter  with  the  double  gate  may  be  ob~ 

SMned  in  similar  fashion.  The  particular  form  chosen  in  writing  U.j 

2 2 

terms  in  the  expressions  for  e_  and  * given  in  the  text  was 

7 * 


the  one  found  moat  oonvonient  for  the  purpoae  of  computation* 

In  the  oaae  of  the  model,  (15)  (for  the  square  law  deteotor) 
and  (16)  (for  the  linear  deteotor  * p » 1)  muat  be  aubstltuted  in 
(8)  and  (9)  to  obtain  d and  d2  , reapeotlvely.  The  expreaalona 
In  the  lntegranda  are  aeeond  degree  polynomials  and  hence  the  Inte- 
grations present  no  dlfflcultlea , However,  several  cases  must  be 
distinguished  with  the  model,  depending  on  the  values  of  06  and  P T, 
since  the  autocorrelation  function  minus  Its  d.c*  value  la  zero  for 
P | x | » 1 • For  the  single  gate  we  muat  consldor  P&  = 1 and 
p 6 > 1,  In  performing  the  Integrations  for  the  double  gate  we 
must  consider  the  following  six  cases i 

1)  P b < 1 

P (T  + 5 ) <1 

2)  pa  < 1 

pT<  1 i p(T  + 6 ) 

3)  pa  < i 

p(T  - a)  <1  = pt 

1*.)  p a < i 

p ( t - a)  = i 

5)  pa  l l 

p(T  - a)  <i  i pt 

6)  pa  £ i 

p(T  - a)  2 i 

2 2 

In  order  to  obtain  by  and  o f In  the  limit  of  Pt  or  p T 

very  large  or  vory  small.  It  is  useful  to  note  that  for  small  x. 
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Q(x)  /— 

0»(x)/^|  from  ( A— 3 ) (A-ll) 

and 

Erf(x)*>-'fr  from  (A-10) 

% 

For  large  x wo  have  (of.  Ref.  A,  p,  222,  Eq.  (8,8)) 

1C 

G'(x)  <— *—  x - (log  x + 1 + C) 

G’(x)^|  (A-12) 

where  C * Euler’s  constant  = 0.5772  . . . 
and  Erf(x)  ^1  . 


As  has  been  noted  earlier,  comparison  may  be  made  with  the  re- 
sults of  Rice  ( Rof , A)  for  noise  alone  in  the  case  of  the  ideal  fil- 
ter, single  and  double  gate  and  in  the  case  of  the  Gaussian  filter., 
single  gato.  For  the  ideal  filter,  single  gate,  the  result  is  given 
in  Ref,  A,  p.  223,  Eq.  (8.10).  Noting  that  where  Rice  uses  <*T  , 

» fb  " fa  and  x^F(x)  *e  use  dy  , W,  &,  P and  G(x),  re- 
spectively, we  may  sot  p =*  0 in  (17),  which  is  then  exactly  four 
times  the  value  given  by  Rice,  as  it  should  be.  (See  comment  in 
paragraph  boginning  on  line  5,  page  13.)  For  the  ideal  filter  with 
double  gato  the  results  in  Ref.  A are  given  on  p.  22lf,  Eq.  (7.3).  in 
the  form  of  an  infinite  aeries  which  may  be  expressed  in  terns  of 
F(x),  as  noted  there,..  Using  Eq.  (6,7)  of  Ref,  A and  noting  that 
••  e>^T  , S,T  of  Ref.  A corrospond  to  d 8 , T,  & , respectively 

of  this  paper,  one  again  finds  that  (after  setting  n®  0 in  (28)), 

2 „ 2 

0 z is  four  times  ° as  given  in  Ref#  A# 

Finally,  the  result  for  the  Gaussian  filtor,  single  gate,  nolee 
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alone  la  given  In  Ref*  A (p.  22$)  in  the  form  of  an  Infinite  series* 


This  aeries  may  be  sunned  aa  shown  below  and  the  result  compared 
with  (16)  for  the  partioular  ease  p « 0, 

From  p*  22$  of  Ref.  A we  have 


dip  * w 


sL  -M 


n ^2n 


0 2 * n*o  (n+1) J(2n+1) 


whore  a * 05 


W 

wo  = 0 


In  the  notation  we  havo  been  using*  Thus,  If  we  define 

«»  - 44  - £ tiiL 

m2 K 2 nsn  t nil ^ • 

thtii 


,2.2n 


W262  n*o  (n+l)l(2n+l) 
co 


• v , _ >n  2n 

[af (a)]  a ^ —ft f .. 

n=o  (n+1). 


or 


a2[af(a>r  --iifSiC.l-i. 

n=o  (n+1); 


n=o 


» 1 - e 


-a2 


I r.  to  grating,  wo  havo 


af(a) 


dt 


J(l-8*‘2)d(-  i] 


V w 

• &A  * 2 1 •-  - 

J O *1 


y/k  Erfta)  + - 


-»2-l 


o:' 


f(a) 


-a2 


Erf  (a)  + Si. 


Substituting  a * 06  y/2ic  , we  have 
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S 


I 


[2  06  y/Tx  Erf(0  6 y/Tx)  ♦ 


e-2  ti02  6 2 


2 

which  is  one  fourth  of  6y  aa  given  in  (18)  with  p » 0,  as  it 
should  bo. 


\ 

i 
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TABLE  II 

Values  of  G(x)  and  CM(x)* 


